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1.  Introduction 


This  report  discusses  the  use  of  infinity  norms  to  solve  the  problem  of  finding  the  minimum 
covering  radius  for  a  set  of  points.  The  minimum  covering  radius  can  be  used  as  a  measure  of 
the  accuracy  of  a  collection  of  shots  or  a  description  of  spall  fragments.  The  infinity  norm 
differs  from  the  1-norm  and  the  2-norm  in  that  it  is  determined  by  a  single  element  of  a  set.  The 
infinity  norm  is  the  largest  deviation  from  a  reference  function  or  point,  the  1-norm  is  the 
average  deviation  (always  positive)  from  the  reference  function,  and  the  2-norm  is  the  average  of 
the  square  root  of  the  square  of  the  difference  between  the  set  of  interest  and  a  reference  set.  The 
2-norm  has  been  extensively  used  as  it  is  commensurate  with  mathematical  analysis.  The 
1-norm  has  become  more  popular  as  computing  has  made  iterative  methods  more  accessible. 

An  appointment  of  infinity  norms  is  that  they  allow  nonprobabilistic  statements  to  be  made. 
These  statements  are  typically  related  to  a  performance  metric  of  some  type  and  have  the 
following  style:  If  the  parameters  are  in  this  range,  then  the  performance  is  better  than  this 
value.  Notice  that  this  form  of  statement  is  typical  of  mathematics  and  is  the  type  of  statement 
required  for  many  artificial  intelligence  database  applications;  that  is,  it  fits  the  framework  of  a 
predicate  calculus.  In  this  report,  the  infinity  norm  is  used  to  make  descriptive  statements  of  the 
form:  All  the  data  is  within  a  circle  of  radius  r,  or  equivalently,  if  the  data  is  from  set  A,  then  it 
is  all  within  a  circle  with  a  radius  of  r.  The  method  developed  can  be  used  for  any  n-dimensional 
space  for  any  performance  metric. 


2.  Method 


A  method  to  find  the  minimum  covering  radius  for  a  set  of  points  was  developed  within  the 
framework  of  infinity  norms.  In  this  case,  the  set  of  controllable  parameters  is  the  center  of  a 
circle.  The  task  is  to  adjust  the  center  in  a  manner  as  to  reduce  the  radius  of  the  circle.  The 
infinity  norm  is  the  greatest  distance  between  each  point  and  the  center.  A  fine  way  to  think  of 
the  algorithm  is  as  an  adjustable  center  removal  process.  For  an  infinity  nonn,  typically  only  one 
point  in  the  data  set  is  pertinent  for  a  given  set  of  conditions.  On  a  given  iteration  only  the  point 
associated  with  the  infinity  norm  and  the  set  of  parameters  are  pertinent  to  reducing  the 
magnitude  of  the  norm.  The  gradient  of  the  distance  function  for  a  particular  point  is  toward  the 
given  point.  For  the  minimum  radius  coverage  problem,  this  amounts  to  moving  the  center 
towards  the  point  associated  with  the  infinity  norm.  The  question  then  is  how  to  determine  the 
magnitude  of  this  displacement.  The  problem  is  simple;  it  is  not  desirable  to  move  so  far  that 
another  point  defines  the  infinity  nonn  and  has  a  larger  magnitude. 
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As  a  first  case,  consider  two  points  and  assume  the  center  is  collocated  at  one  of  the  points. 
Assuming  distance  is  being  used  as  the  metric;  one  point  will  have  a  norm  equal  to  the  distance 
between  the  two  and  the  other  will  have  a  norm  of  zero.  If  the  magnitude  of  the  difference  in  the 
norms  is  used  as  an  adjustment  for  the  center,  the  process  will  oscillate  between  the  two  points 
and  never  converge.  This  oscillation  occurs  along  the  direction  connecting  the  two  points,  and  is 
eliminated  by  using  0.5  as  a  multiplier  of  the  magnitude.  In  this  situation,  when  the  midpoint 
between  the  two  points  is  reached,  there  are  two  points  with  the  same  norm.  Any  change  to  the 
circles  center  from  the  midpoint  can  be  thought  of  as  a  component  along  the  connecting  segment 
and  a  component  orthogonal  to  this  segment.  Any  change  along  the  connecting  segment  will 
increase  the  infinity  norm  as  one  distance  (metric)  must  increase.  Changes  in  the  orthogonal 
component  will  cause  both  distances  (metrics)  to  increase.  For  distances,  the  Pythagorean 
Theorem  guarantees  this  increase;  for  other  metrics,  the  Cauchy  Schwartz  inequality  can  be 
invoked.  So  a  deviation  from  the  midpoint  along  either  component  predicts  an  increase  of  the 
infinity  norm.  This  suggests  an  exit  condition  for  an  algorithm  can  be:  If  the  center  point  is  at 
the  midpoint  of  two  points  and  the  norms  associated  with  all  the  other  points  are  less  than  these 
two  norms,  then  the  search  is  done. 

Next,  consider  an  exit  condition  for  a  set  of  three  points,  if  the  third  point  is  within  the  circle 
defined  by  the  two  most  distant  points  the  situation  reverts  to  the  previous  case;  however,  the 
third  point  can  keep  the  center  from  reaching  the  midpoint  of  the  two  most  distant  points.  If  the 
center  is  located  so  that  all  three  points  are  on  the  circumference  of  the  same  circle  then  the 
algorithm  should  exit  as  there  is  no  basis  for  improvement  in  reducing  the  infinity  norm. 
Consider  the  situation  where  the  algorithm  has  found  two  points  that  have  the  same  infinity 
nonn.  At  this  point  for  the  distance  metric,  the  center  of  the  circle  must  be  on  the  perpendicular 
bisector  of  the  two  points.  To  decrease  the  norms  associated  with  these  two  points,  the  center 
must  move  toward  their  midpoint;  during  this  process,  it  is  possible  for  a  third  point  to  prevent 
further  decrease  when  it  obtains  the  same  value  as  the  other  two  points.  In  this  situation,  the 
center  will  be  at  the  intersection  of  the  perpendicular  bisectors. 

Based  on  the  former  discussion,  an  algorithm  to  find  the  minimum  covering  radius  was 
developed.  The  basic  step  is  to  move  towards  the  point  with  the  largest  distance  from  the  center. 
This  process  is  repeated  until  the  difference  between  the  two  largest  distances  is  within 
tolerances  defined  by  the  difference  in  the  distances  between  the  second  and  third  largest 
distances.  If  the  criteria  is  met,  the  process  moves  toward  the  midpoint  of  the  segment  defined 
by  the  two  points  with  the  largest  distances.  If  the  midpoint  is  reached,  the  algorithm  terminates; 
if  three  points  have  equal  distances  (the  same  norm),  then  the  process  terminates.  The  algorithm 
uses  a  sort  routine  to  order  the  distance  metric  and  then  a  routine  to  find  the  index  of  a  point  with 
a  given  value.  Both  of  these  functions  are  available  as  MATLAB  commands. 
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3.  Implementation 


The  ideas  were  implemented  through  a  graphical  user  interface  (GUI)  in  MATLAB  shown  in 
figure  1.  There  are  three  buttons  on  the  GUI.  The  first,  called  “Get  Data,”  queries  the  user  for  a 
data  file.  The  file  is  assumed  to  contain  two  rows  of  entries  for  the  X  and  Y  values  of  the  data. 

If  the  data  is  stored  from  MATLAB,  it  must  be  in  the  first  element  using  the  save  command.  For 
example,  if  the  data  is  in  a  variable  named  dat,  then  the  command,  “save  info  dat,”  will  place  the 
information  in  the  variable  dat  into  the  file,  “info.mat.”  After  the  file  is  opened,  the  GUI  will 
display  a  graph  of  the  data.  To  continue,  the  user  must  specify  the  max  number  of  replications  to 
be  made.  There  is  a  slide  bar  to  adjust  this  value;  alternatively,  the  user  can  enter  a  number  in 
the  text  box  above  the  slide  bar.  After  the  maximum  number  of  replications  is  set,  the  user  then 
clicks  the  “Adjust  Center”  button.  The  ensuing  graph  shows  the  movement  of  the  center.  The 
number  of  iterations  to  meet  termination  requirements  is  displayed  in  the  replications  text  box. 
The  user  can  click  the  “Final  Graph”  button  to  see  the  end  product.  Displayed  on  the  graph  will 
be  the  data  points  in  blue,  the  center  of  the  circle  as  a  red  plus,  and  the  circle  in  green.  On  the 
GUI  under  the  title  “Center,”  the  center  of  the  circle  will  be  displayed;  under  the  title  “Radius,” 
will  be  the  minimum  covering  radius.  Finally,  if  the  user  wants  a  figure  for  reports  without  the 
GUI  interface,  the  menu  item  “AxestoFigure”  can  be  clicked.  This  will  open  a  figure  window 
with  the  final  graph  using  the  file  name  as  the  title  and  the  center  location  and  radius  printed 
below  the  x-axis.  Properties  can  be  adjusted  using  the  figure  menu  item  axes  properties  (under 
the  edit  menu).  After  this,  the  figure  can  be  copied  and  moved  to  any  document.  The  MATLAB 
program  listing  is  included  in  the  appendix. 


Figure  1.  Graphics  user  interface  for 
covering  circle  algorithm. 
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The  starting  point  for  the  center  was  set  to  the  mean  of  the  locations.  This  worked  for  all  the 
data  sets  investigated.  The  pseudo-code  is  contained  in  figures  2  and  3. 


Main  Function 

Data  is  an  array  of  points 

Choose  Center  as  Average  of  Data 

lastCenter=Center+2*tolerance 

For  i=l  till  #  of  iterations 

If  distance  between  consecutive  centers  <  tolerance 
Done 

Else 

Call  MoveCenter 

Endif 

Endfor 

Done 


Figure  2.  Main  function  pseudo-code. 


Function  MoveCenter 

//comment  costfunction  evaluates  the  data  based  on  the  metric 
//comment  sort  puts  the  data  in  descending  order  high  to  low 
N  e  wdata=data-center 
D2=costfunction(N  ewdata ) 

Ds=sort(D2) 

If  Ds(  1  )=Ds(2)=Ds(3) 

Done 

Endif 


Figure  3.  Function  MoveCenter  pseudo-code. 


It  is  prudent  to  mention  the  following  situation  as  an  example  of  the  algorithm  converging  to  an 
improper  solution.  Consider  three  points  on  a  circle  within  a  small  arc.  If  the  algorithm  selects 
the  center  of  this  circle  on  any  iteration  it  will  terminate.  This  is  an  undesirable  situation  since 
the  distance  between  the  points  will  all  be  less  than  the  radius  of  this  circle;  this  is  illustrated  in 
figure  4.  Although  this  situation  did  not  occur  using  the  mean  of  the  locations  as  the  starting 
point  for  the  center,  it  did  occur  in  one  testing  situation.  Observation  of  the  graph  of  the  solution 
can  be  used  to  rule  out  this  case.  To  test  a  solution  it  would  be  possible  to  restart  the  algorithm  at 
the  average  of  the  three  points  defining  the  final  circle  and  then  determine  if  the  algorithm 
converges  to  the  same  solution. 
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Figure  4.  An  inappropriate  covering  circle. 


4.  Conclusion 


It  is  certainly  possible  to  develop  a  more  complex  algorithm  utilizing  more  detailed  information. 
For  example,  the  distant  matrix  could  be  calculated  for  the  set  of  n  points;  this  will  take  0.5n2  -  n 
distance  calculations.  The  Elzinga-Heam  method1  is  a  geometric  algorithm  that  solves  this 
problem.  The  method  developed  is  conceptually  straightforward  as  it  is  based  on  the  definition  of 
the  infinity  norm.  Unlike  methods  associated  with  linear  programming,  this  method  does  not 
require  linear  metrics.  The  algorithm  worked  well  for  the  data  sets  investigated,  sometimes 
converging  in  three  iterations;  however,  in  some  cases,  there  were  hundreds  of  iterations.  For 
specific  metrics,  it  would  be  possible  to  use  directional  derivatives  to  improve  the  convergence 
of  the  process.  The  overall  design  is  based  on  defining  an  improvement  step  to  be  repeated  until 
the  state  of  the  process  fulfills  a  specific  criterion.  Infinity  norms  offer  a  theoretic  framework  for 
algorithm  development. 


1 


Elzinga,  D.  J.;  Hearn,  D.  W.  Geometrical  Solutions  for  Some  Minimax  Location  Problems.  Transportation  Science  1972,  6, 


379-394. 
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Intentionally  left  blank. 
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Appendix.  MATLAB  Program  Listing 


This  appendix  appears  in  its  original  form,  without  editorial  change. 
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function  varargout  =  ccgui  (varargi  n) 

%  CCGUI  M-  f  i  I  e  for  ccgui  .  f  i  g 

%  CCGUI,  by  itself,  creates  a  new  CCGUI  or  raises  the  existing 

%  singleton*. 

% 

%  H  =  CCGUI  returns  the  handle  to  a  new  CCGUI  or  the  handle  to 

%  the  existing  singleton*. 

% 

%  This  function  finds  the  radius  of  a  the  covering  circle  of  a  set 
%  of  points  ( x ,  y ) 

%  The  button  get  data  is  used  to  open  a  file  containing  the  data 

%  The  button  adjust  data  is  used  in  conjunction  with  the  number  of 

%  iterations  to  cycle  in  an  attempt  to  find  the  center  of  the  circle 

%  with  the  smallest  covering  radius 

% 

% 

%  CCGUI  (' CALLBACK' ,  hObj  ect ,  event  Dat  a,  handl  es,  ...  )  calls  the  local 

%  function  named  CALLBACK  in  CCGUI  ,  M  with  the  given  input  arguments, 

% 

%  CCGUI  ( '  Property'  , 1  Val  ue'  ,  .  .  .  )  creates  a  new  CCGUI  or  raises  the 
%  existing  singleton*.  Starting  from  the  left,  property  value  pairs  are 

%  applied  to  the  GUI  before  ccgui  Openi  ngFcn  gets  called.  An 
%  unrecognized  property  name  or  invalid  value  makes  property  application 

%  stop.  All  inputs  are  passed  to  ccgui  Openi  ngFcn  via  varargin. 

%  " 

%  *See  GUI  Options  on  GUIDE'S  Tools  menu.  Choose  "GUI  allows  only  one 
%  instance  to  run  (singleton)". 

% 

%  See  also:  GUIDE,  GUI  DATA,  GUI  HANDLES 


%  Edit  the  above  text  to  modify  the  response  to  help  ccgui 


%  Last  Modified  by  GUIDE  v  2 .  5  04-  Feb-  2008  1  1:3  2:0  2 
%A  n  d  r  e  w  T  h  o  mp  s  o  n 
%2  0  0  8  US ARL 


gui 
g  u  i 


Begin  initialization  code  -  DO  NOT  EDIT 


Singleton  =  1; 
'State  =  struct! 


i  f 
end 


gui  _ Name1  , 
gui  _Si  ngl  e t  o n 1 
1  gui  _ O p e n i  ngFcn 
'  gui  _Out  put  Fen' 

'  gui  _Layout  Fen' 

1  gui  _  C  a  I  I  back1  , 
nar  gi  n  &&  i  s c  ha r  ( var  ar  gi  n { 1 }) 
gui  _State.  gui  _Cal  I  back  =  s  t  r  2f  u  nc  ( v  a  r  a  r  g  i  n  { 1 } ) 


mf  i 
g  u  i  _  S 
@c  c  g  u 
@c  c  g  u 
M  , 
[]); 


e  n  a  me ,  ... 
ngl  et  on, 
Openi  ngFcn, 
"Out  put  Fen, 


if  n  a  r  g  o  ut 

[varargout{l:  nargout}]  =  gui  ma  i  n  f  c  n  ( g  u  i  State,  varargi  n{:  }); 

else 

gui_mainfcn(gui  State,  varargi  n{:  }) ; 
end 

%  End  initialization  code  -  DO  NOT  EDIT 


%  -  -  -  Executes  just  before  ccgui  is  made  visible. 

function  ccgui  Openi  n  g  F  c  n  (  hObj  ect ,  eventdata,  handles,  varargin) 

%  This  function  has  no  output  args,  see  Out  put  Fen. 

%  hObj  ec  t  handle  to  figure 

%  eventdata  reserved  ■  to  be  defined  in  a  future  version  of  MATLAB 
%  handles  structure  with  handles  and  user  data  (see  GUI  DATA) 

%  varargin  command  line  arguments  to  ccgui  (see  VARARGIN) 

%  C  h  o  o  s  e  default  c  o  mma  n  d  line  output  for  ccgui 
ha ndl es . out  put  =  h Ob j  ect ; 

%  U  p  d  a  t  e  handles  structure 
gui  data)  hObj  ect,  handles); 

%  UIWAIT  makes  ccgui  wait  for  user  response  (see  UIRESUME) 

%  ui  wai  t(handl  es.fi  gurel); 


%---  Outputs  from  this  function  are  returned  to  the  command  line, 
function  varargout  =  ccgui  Out  put  Fcn(  hObj  ect ,  eventdata,  handles) 

%  varargout  cell  array  for  returning  output  args  (see  VARARGOUT); 

%  hObj  ec  t  ha  ndl  e  t  o  f  i  gur  e 

%  eventdata  reserved  ■  to  be  defined  in  a  future  version  of  MATLAB 

%  handles  structure  with  handles  and  user  data  (see  GUI  DATA) 
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%  Ge t  default  c o mma nd  line  output  from  handles  structure 
varargout{l}  =  handles,  output; 


%  -  -  -  Executes  on  button  press  in  pushfile. 

function  pushfile  Cal  I  back(hObj  ect,  eventdata,  handles) 

%  hObj  ect  handfe  to  pushfile  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of  MATLAB 
%  handl  es  structure  with  handles  and  user  data  (see  GUI  DATA) 

[fi  I  ename,  path]  =ui  getfi  I  e( 1  *.  *'  ) ; 
if  i  sequal  (fi  I  ename,  0) 
d i  s p ( 1  no  change1  ) 

else 

set  (  handl  es.textpath,  1  stri  ng1  ,  path); 
setj  handl  es.textfi  I  e,  1  stri  ng1  ,  f  i  I  ename); 
d  a  t  =1  o  a  d  ( f  u  I  I  fi  I  e  ( pat  h ,  fi  I  ename)); 

if  i  s  s  t  r  u  c  t  ( d  a  t )  %i  f  saved  as  a  *  ma  t  it  is  loaded  as  a  structure  take  the  first 
field  as  data 

fnames=fi  el  dnames(dat); 
dat  =d a t .  ( f  na  mes { 1 } ) ; 
end 

[  m,  n ]  =s  i  z  e(  da t ) ; 
i  f  n  ==2 

dat  =d  a  t 1 ; 
end 

axes( handi  es.  axesl) ; 
hold  off 

pi  ot  ( dat  (  1,  :  ) ,  dat  ( 2,  :  ) ,  1  b*'  ) 

set  ( handl  es.  pushf  i  1  e, '  User  Dat  a 1 , dat) ; 

%s  e  t  ( handl  e  s .  t  e  x  t  C  e  n  t  e  r ,  1  stri  n  g 1  ,  1  1  ) ; 

%s  e  t  ( handl  es.textRadi  us,  1  stri  ng1  ,  1  1  j ; 
end 


function  edi  1 2  _  C  a  I  I  back)  hObj  ect ,  eventdata,  handles) 

%  hObj  ect  handle  to  e  d  i  1 2  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of  MATLAB 
%  handl  es  structure  with  handles  and  user  data  (see  GUI  DATA) 

%  Hints:  get  ( hObj  ect ,  'String')  returns  contents  of  e  d  i  1 2  as  text 
%  str2double(get(h  Object, 'String'))  returns  contents  of  e  d  i  1 2  as  a  double 


%---  Executes  during  object  creation,  after  setting  all  properties, 
function  edi  t2_CreateFcn(  hObj  ect,  eventdata,  handles) 

%  hObj  ect  handle  to  e  d  i  1 2  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of  MATLAB 
%  handles  empty  -  handles  not  created  until  after  all  CreateFcns  called 


%  Hint:  edit  controls  usually  have  a  white  background  on  Windows, 
%  See  I  SPC  and  COMPUTER. 

if  ispc  &&  i  sequal  (getfhObi  ect, '  BackgroundCol  or'  ), 
get  (  0,  '  defaul  tUi  control  BackgroundCol  or'  )) 
set  (  hObj  ect,  '  BackgroundCol  or'  ,  '  whi  t  e '  ); 
end 


%  -  -  -  Executes  on  slider  move  me  nt. 

function  si  i  d  e  r  2  Cal  I  b  a  c  k  ( hObj  ect ,  eventdata,  handles) 

%  hObject  handle  to  si  i  d  e  r  2  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of  MATLAB 
%  handl  es  structure  with  handles  and  user  data  (see  GUI  DATA) 

%  Hints:  get  ( h  Ob  j  ect ,  1  Val  u  e '  )  returns  position  of  slider 

%  get  ( hObj  ect ,  '  Mi  n' )  and  get  ( hObj  ect ,' Max' )  to  determine  range  of  slider 

set  ( handl  es.  edi  1 2 , '  stri  ng'  ,  n  u  m2  s  t  r  ( f  I  oor( get (  handl  es.  si  i  d  e  r  2 ,  '  Val  u  e ' )))); 

%■■■  Executes  during  object  creation,  after  setting  all  properties, 
function  si  i  d  e  r  2  Cr  e  at  e  F  c  n  ( hOb  j  ec  t ,  eventdata,  handles) 

%  hObj  ect  handle  to  si  i  d  e  r  2  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of  MATLAB 
%  handles  empty  -  handles  not  created  until  after  all  CreateFcns  called 

%  Hint:  slider  controls  usually  have  a  light  gray  background. 

i  f  i  sequal  (get(  hObj  ect,  '  BackgroundCol  or'  ) ,  get(  0,  '  defaul  tUi  control  BackgroundCol  or1  )  ) 
set(  hObj  ect,  '  BackgroundCol  or'  ,  [  .  9  .9  .9]); 
end 
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%---  Executes  on  button  press  in  pbAdj  ustCi  rcl  e. 

function  pbAdj  ustCi  rcl  e_Cal  I  back(  hObj  ect,  eventdata,  handles) 

%  h  Ob  j  ect  handle  to  pbAdj  us  t  Ci  r  c  I  e  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of  MAT L A B 
%  handles  structure  with  handles  and  user  data  (see  GUI  DATA) 

da t  =get  ( handl  es.  pushfi  I  e,  '  User  Da  t  a 1  ) ; 
n=str2num(get(  handl  es.  edi  1 2 ,  1  stri  ng1  )); 
c  =[  me  a  n  ( d  a  t  ( 1 ,  :  ));  me  a  n  ( d  a  t  ( 2 ,  :  ))  ]; 
cl  i  s t  =c ; 
for  i  =1:  n 
cl  ast  =c; 

[  c,  r  ]  =moveCent  er  ( dat ,  c) ; 
cdi  f  =c- cl  ast ; 
ncdi  f  =nor  m(  cdi  f ) ; 
if  ncdi  f  <.  00000001 
break; 
end 

cl  i  st  =[  cl  i  st ,  c] ; 
end 

set  ( handl  es.  edi  1 2 ,  1  stri  ng1  ,  n  u  m2  s  t  r  ( i  )); 
a x e s ( handl  es.  axesl) ; 

pi  ot  ( cl  i  st  (  1,  ;  ) ,  cl  i  st  ( 2, :  ) ,  1  b- '  ,  cl  i  st  (  1, :  ) ,  cl  i  st  ( 2, ;  ) ,  '  r .  1  ,  c( 1) ,  c(  2) ,  1  r +'  ) 
set(handl  es.  textCenter,  1  stri  ng1  ,  [num2str(c(l),  4),  1  ,  ',num2str(c(2),4)]); 
set  ( handl  es.  t  ext  Radi  us,  '  stri  ng'  ,  n  u  m2  s  t  r  ( r ,  5 ) ) ; 
c  r  =[  c ;  r ] ; 

set  ( handl  es.  pbAdj  ustCi  rcl  e,  '  UserData1  ,  cr); 


%  -  -  -  Executes  on  button  press  in  pushbutton3. 

function  pushbutton3  Cal  I  back(  hObj  ect,  eventdata,  handles) 

%  hObj  ect  handle  to  pushbutton3  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of  MAT  LAB 
%  handl  es  structure  with  handles  and  user  data  (see  GUI  DATA) 

dat  =g  e  t  (  handl  es.  pushfi  I  e,  1  UserData1  ) ; 
cr=get(  handl  es.  pbAdj  ustCi  rcl  e,  1  UserData1  ); 
t  het  a=0:  .  01;  2 * p i  +.  01; 
a  x  e  s  (  handl  es.  axesl) ; 

pi  ot(cr(3)*cos(theta)+cr(l),  c  r  ( 3 )  *  s  i  n(theta)+cr(2),  '  g'  ) 
hoi  d  on 

pi  ot ( dat ( 1,  :  ) ,  dat  ( 2,  : ) ,  1  b*1  ) 
pi  ot  ( c r (  1) ,  c r  ( 2) ,  1  r  +'  ) 


% . 

function  Untitled  1  Cal  I  back(  hObj  ect,  eventdata,  handles) 

%  hObj  ect  handTe’to  Unt  i  1 1  ed_l  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of  MAT  LAB 
%  handl  es  structure  with  handles  and  user  data  (see  GUI  DATA) 


% . 

function  Untitled  3  Cal  I  bac  k(  hObj  ect ,  eventdata,  handles) 

%  hObj  ect  handTe'to  Untitled_3  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of  MATLAB 
%  handles  structure  with  handles  and  user  data  (see  GUI  DATA) 


% 

function  Untitled  4  Ca  I  I  ba  c  k  ( h  Obj  ec  t ,  eventdata,  handles) 

%  hObj  ect  handTe'to  Untitled_4  (see  GCBO) 

%  eventdata  reserved  ■  to  be  defined  in  a  future  version  of  MATLAB 
%  handles  structure  with  handles  and  user  data  (see  GUI  DATA) 


% . 

function  Untitled  5  Cal  I  back(  hObj  ect,  eventdata,  handles) 

%  hObj  ect  handTe’to  Untitled_5  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of  MATLAB 
%  handl  es  structure  with  handles  and  user  data  (see  GUI  DATA) 


% . 

function  Untitled  2  Cal  I  bac  k(  hObj  ect ,  eventdata,  handles) 

%  hObj  ect  handTe'to  Untitled_2  (see  GCBO) 

%  eventdata  reserved  -  to  be  defined  in  a  future  version  of  MATLAB 
%  handles  structure  with  handles  and  user  data  (see  GUI  DATA) 


% 

function  ui  toggl  etool  3_CI  i  ckedCal  I  back(  hObj  ect,  eventdata,  handles) 
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%  h  Ob  j  ect  handle  to  ui  toggl  etool  3  (see  GCBO) 

%  eventdata  reserved  ■  to  be  defined  in  a  future  version  of  MAT  LAB 
%  handles  structure  with  handles  and  user  data  (see  GUI  DATA) 


% . 

function  File  menu  Cal  I  b  a  c  k  ( hObj  ect ,  eventdata,  handles) 

%  h  0  b  j  e  c  t  handle  to  File  menu  (see  GCBO) 

%  eventdata  reserved  -  to  Be  defined  in  a  future  version  of  MATLAB 

%  handl  es  structure  with  handles  and  user  data  (see  GUI  DATA) 

newf igurehandle  =  figure 

copyobj  ( handl  es.  axes  1 ,  newfi  gurehandl  e) 

t=get(  handl  es.  textfi  I  e,  1  stri  ng'  ) ; 

t  =[ '  F i  I  e :  1  ,  t  ] ; 

c=get(  handl  es.  textCenter,  1  stri  ng1  ); 
xa  =[ 1  Center:  ( '  ,  c,  '  ) 1  ] ; 
r=get(  handl  es,  t  ext  Radi  us,  1  stri  ng1  ); 
x a  =[  xa,  1  Radi  us:  1  ,  r] ; 
t  i  1 1  e  ( t ) 
xl  abel ( xa) 

pri  nt(newfi  gurehandl  e,  '  -  d  me  t  a 1  ) 

%d  e  I  e  t  e  ( newfi  gurehandl  e ) 


function  [center,  r]  =moveCenter(  dat,  center) 

%d  a  t  is  a  2  d  i  me  n  s  i  o  n  a  I  set  of  data 
%c  is  the  a  s  s  u  me  d  center 

% 

%t  h  e  algorithm  works  by  moving  the  center  toward  the 
%d  a  t  a  point  with  the  greatest  distance  it  is  done  when 
%two  or  three  points  are  all  at  the  same  approximate  distance 
%f  r  o  m  the  center 

%t  h  i  s  code  is  just  one  step  so  it  needs  to  be  called  by  a  routine  that 
%decides  when  to  stop  adjusting 

%A  n  d  r  e  w  T  h  o  mp  s  o  n 

%b  a  s  i  cal  I  y  an  infinity  norm  minimization  routine 


[  r ,  c] =si  ze( dat ) ; 
if  c  ==2 

dat  =d  a  t ' ; 
n  =r ; 

else 

n=c; 

end 

[  r ,  c ] =s  i  z  e( cent  er ) ; 
if  c  ==2 

center  =cent  er 1 ; 
end 

z=repmat  ( cent  er,  1,  n); 
newdat  =d a t  -  z ; 
ndat  2  =n e wd a t .  A2; 
d 2  =s u m(  ndat  2) ; 
d  2  =s  q  r  t  ( d  2 ) : 

ds=sort(d2,  2,  1  descend1  ) ; 

delta  =.  5*(ds(l)-ds(2));  %.  5  minimizes  the  convergence  oscillation 

i  =f  i  nd(  d 2  =  =  ma x (  d 2 )  j ; 

%more  than  two  points  on  circumpherence  so  done 
if  I  engt  h( i ) >2 
r  =ds ( 1) ; 
return 
end 

%i  f  I  engt  h(  i )  ==2 

if  (  (  ds(  1) - ds( 2) )  <  .  2*(  ds( 2) - ds( 3) )  ) 
if  length(i)  <2 

i=[i  f  i  nd(  d2==ds(  2) )  ] ; 
end 

mi  d=(  dat  (  :  ,  i  ( 1) )  +d a t  ( :  ,  i  ( 2) ) )  /  2; 
if  center  ==mi  d 

%o  u  t  = '  centered' 
r  =ds(  1) ; 
return: 
end 

dc  1  =mi  d  -  c  en t  e r ; 

%dc  =dc  1/  nor  m(  dc  1) ; 
del  t  a  =ds ( 1)  -  ds ( 3) ; 


11 


if  del  t  a  >n  o  r  m(  dc  1) 
del  t  a  =1 ; 
end 

cent  er  =cent  er  +d  e I  t  a  *  d  c  1 ; 

else 

dc  =d a t  ( :  ,  i  )  •  c ent  er ; 
dc  =d c /  nor  m(  dc ) ; 
cent  er  =cent  er +  d  e I  t  a  *  d  c ; 
end 

z  =r  epmat  ( cent  er ,  1,  n) ; 
newdat  =d  a  t - z ; 
ndat  2  =n e wd a t .  *2; 
d2=sum(  ndat  2) ; 
d  2  =s  q  r  t  (  d  2 ) ; 

d  s  =s  o  r  t  ( d  2 ,  2,  'descend'); 
r  =ds(  1) ; 


12 


NO.  OF 

COPIES  ORGANIZATION 


1  DEFENSE  TECHNICAL 
(PDF  INFORMATION  CTR 
ONLY)  DTICOCA 

8725  JOHN  J  KINGMAN  RD 
STE  0944 

FORT  BELVOIR  VA  22060-6218 

1  US  ARMY  RSRCH  DEV  & 
ENGRG  CMD 
SYSTEMS  OF  SYSTEMS 
INTEGRATION 
AMSRD  SS  T 
6000  6TH  ST  STE  100 
FORT  BELVOIR  VA  22060-5608 

1  DIRECTOR 

US  ARMY  RESEARCH  LAB 
IMNE  ALC  IMS 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1197 

1  DIRECTOR 

US  ARMY  RESEARCH  LAB 
AMSRD  ARL  Cl  OK  TL 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1197 

1  DIRECTOR 

US  ARMY  RESEARCH  LAB 
AMSRD  ARL  Cl  OK  T 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1197 


ABERDEEN  PROVING  GROUND 


1  DIR  USARL 

AMSRD  ARL  Cl  OK  TP  (BLDG  4600) 


13 


NO.  OF  NO.  OF 

COPIES  ORGANIZATION  COPIES 

1  DREXEL  UNIV 

DEPT  OF  MECHL  ENGRG 
B  CHANG 
3141  CHESTNUT  ST 
PHILADELPHIA  PA  19104 

1  DIRECTOR 

US  ARMY  RSRCH  LAB 
AMSRD  ARL  Cl  NT 
R  PRESSLEY 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1197 


ABERDEEN  PROVING  GROUND 

2  COMMANDER 
USAATC 
TEDT  AT  ADR 
B  G1LLICH 
S CLARK 

400  COLLERAN  RD 
TRAILER  T1 
APG  MD  21005-5059 

1  COMMANDER 
USAATC 
TEDT  AT  ADF 
S  NOVAK 
400  COLLERAN  RD 
APG  MD  21005-5059 

21  DIR  US  ARL 

AMSRD  AAR  AEF  T 
M  ANDRIOLO 
AMSRD  ARL  Cl  CT 
B  BODT 
RKASTE 

AMSRD  ARL  SL  BD 
J  COLLINS 
L  MOSS 

AMSRD  ARL  SL  BW 
P  GILLICH 
AMSRD  ARL  WM  BA 
R  MCGEE 
T  BROWN 
T  HARKINS 
MILG 

AMSRD  ARL  WM  BF 
J  WALD 
J  WALL 
D  WEBB 


ORGANIZATION 

M  ARTHUR 
A  THOMPSON  (4  CPS) 
B  FLANDERS 
R  PEARSON 
B  OBERLE 


14 


